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Abstract 

Context. Markarian 509 and Markarian 841 are two bright Seyfert 1 galaxies with X-ray spectra characterised by a strong soft excess 
and a variable Fe Ka line, as shown by several X-ray observatories in the past. 

Aims. We report an analysis and modelling of new Suzaku observations of these sources, taken between April and November, 2006, 
for Mrk 509, and between January and July, 2007, for Mrk 841, for a total exposure time of » 100 kiloseconds each. Data from 
XIS and HXD/PIN instruments, going from 0.5 to 60 keV, represent the highest spectral resolution simultaneous broad-band X-ray 
spectrum for these objects, and provide the strongest constraints yet on the origin of the soft excess emission. 

Methods. We fitted the broad-band spectrum of both sources with a double Comptonisation model, adding neutral reflection from 
distant material and a two-phase warm absorber. We then studied the two competing models developed to explain the soft excess in 
terms of atomic processes: a blurred ionised disc reflection and an ionised absorption by a high velocity material. 
Results. When fitting the data in the 3-10 keV range with a power law spectrum, and extrapolating this result to low energies, a soft 
excess is clearly observed below 2 keV, although its strength is weak compared to previous observations of both sources. A moderate 
hard excess is seen at energies higher than 10 keV, together with a neutral Fe Ka narrow emission line at E » 6.4 keV and a broad Fe 
emission line. For Mrk 509, the broad Fe emission line is required in all the three physical models to ensure a good fit to the data: this 
finding suggests that the blurred reflection model correctly describes the soft excess, but that it underestimates the broad Fe emission 
line. For the smeared absorption model, this suggests instead that the continuum spectrum absorbed by the outflowing gas should 
indeed contain a reflected component. For Mrk 841, all three models that we tested provide a good fit to the data, and we cannot rule 
out any of them. A broad emission line is required in the double Comptonisation and smeared absorption models, while the blurred 
reflection model consistently fits the broad-band spectrum, without adding any extra emission-line component. 
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1. Introduction 

The origin of the soft X-ray excess emission is still nowadays, 



about three decades after the discovery (Pravdo et al. 1981 
|Arnaud et al.| fl985), one of the major open questions in AGN 
(active galactic nuclei) research (e.g. |Turner & Miller 2009). 
This component appears as a featureless excess of emission 
above the low energy extrapolation of the 2-10 keV best-fit 
power law. Historically associated with the high energy tail of 



the accretion-disc black-body radiation, many autho rs (|Czerny 



etalT2 003; 
et al.||2005 



Gierlihski & Done|200 4 hereafter GD(M; [Piconceil7 



Crummy et al.| [2006 hereafter C06) have shown 
that modelling with a thermal continuum infers a characteristic 
"temperature" higher than expectations (Shakura & Sunyaev] 
|1976] l that remains remarkably constant across a range of AGN 
despite a wide spread in black hole mass and AGN luminosity. 
Moreover, in bright and variable AGN the soft excess does 
not follow the expected black-body luminosity-temperature 
relation ( |Ponti et al.||200 6). Various studies have shown that 
the ratio of the soft excess, at 0.5 keV, to the extrapolation of 



the high-energy power-law emission has a very small scatter 
( ]Piconcelli et al.||2005"| |Miniutti et~al"| |2009| l, which is very 
different from what is observed in galactic black holes in the 
bright soft state, dominated by the disc black-body emission 
(e.g. |Done et al"||2007| ). These findings triggered the search for 
a different origin of the soft excess. The first idea investigated 
was the Comptonisation that might occur in the upper layer of 
the accretio n disc ( |Czerny et"al"||2003| GD04; |Sobolewska &| 
Done 2007). This model may explain the featureless shape of 
the soft excess and its high temperature; the feedback between 
the coronal disc irradiation and the disc skin may provide a 
mechanism to stabilize both the soft excess temperature and 
the ratio relative to the power-law emission. On the other hand, 
the observed constancy of the soft excess temperature seems 
to suggest a nature tied to atomic processes. Two possible 
alternatives are reflection and absorption. If the upper layer 
of the accretion disc is ionized, the disc Compton reflection 
component will contain many X-ray lines (Ross & Fabian 2005 ) 
that will be broadened by the relativistic motion of the material 
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Table 1. Suzaku observation log 



Mrk 509 




ID 


Aim Point 


Start 


Exposure (s) 


obs.l 


701093010 


HXD 


2006-04-25 


24580 


obs.2 


701093020 


XIS 


2006-10-14 


25960 


obs.3 


701093030 


XIS 


2006-11-15 


24450 


obs4 


701093040 


XIS 


2006-11-27 


33090 


tot 








108080 


Mrk 841 




ID 


Aim Point 


Start 


Exposure (s) 


obs.l 


701084010 


HXD 


2007-01-22 


51790 


obs.2 


701084020 


HXD 


2007-07-23 


50930 


tot 








102720 



on the surface of the accretion disc, forming a featureless 
continuum (C06). Ionized absorption processes can also imprint 
their main features in the soft X-ray band. If these absorbers are 
moving relativistically (for example, in the form of a disc wind) 
they may, as well, produce a featureless curved continuum and 
reproduce the soft excess. Nevertheless, numerical simulations 
indicate that, to reproduce the excess, the velocities of such 
a wind have to be extremely high (Schurch & Done 2008), 
suggesting that the flow has to be clumpy and/or only partially 
covering the source, or associated with a magnetically driven jet. 



These different models are degenerate in the 0.5-10 keV 
band but still predict different amounts of hard X-ray emission. 
In addition, they predict different variability behaviours in the 
soft and hard energy band and, thus, the measurement of the 
high energy flux may help us to determine the physical origin of 
the soft excess. To this aim, with Suzaku we observed Markarian 
509 (Mrk 509), the brightest Seyfert 1 galaxy in the hard X- 
ray sky (17-60 keV) with an important soft excess but no strong 
warm absorber component ( |Sazonov et al.|20 07), and Markarian 
841 (Mrk 841), one of the first AGN in which a soft excess has 
been discovered ( Arnaud et al. 1985). Information about the four 
Mrk 509 observations and the two observations of Mrk 841 are 
given in Table [TJ 



Located at a redshift of z = 0.0344 ( |Fisher et ah] 



1995), 



Cooke 



Mrk 509 was first observed at X-ray energies by Ariel V 
|et al.|[T978] ), and then studied by several X-rays telescopes 
Using g£AO-7, |Singh et aL] ( |1985[ ) detected its soft excess for 
the first time, a detection that was later confirmed by EXOSAT 
(Mor ini et al.||198"7] l, in addition to the detection of a Fe Ka 
emission line. Using Ginga, |Singh et al.| ( |1990[ ) revealed a 
hardening of the power law at energies greater than 10 keV. 
BeppoSAX observations ( |Perola et al.||2QQ0] > showed a power- 
law cut-off at 70 keV and revealed information about the warm 
absorber (Dadina et al. 



depth by Smith et al. 2007 ) and Detmers et al. 



2005 1 , which since has been studied in 

with 



XMM-Newton/RGS data. Ponti et al. (2009) presented the most 



complete study to date of the Fe K complex, using XMM- 
Newton data and the same Suzaku observations we present here. 
For this reason, we focus on the broad-band spectrum only, and 
no further analysis of emission lines will be done. Recently, 
|Kaastra et al.| ( |20TT] > carried out a multi -wavelength campaign 
of this source, using data from seven different satellites and 



observatories (from infra-red to hard X-rays, see Kaastra et al. 
( |201 l| l and references therein). 

The soft excess of the bright Seyfert 1 galaxy Mrk 841 
(z = 0.0365) was first detected by EXOSAT ( [Arnaudet al.|1985| l. 




soo 

200 



Right ascension 

Figure 1. Suzaku/XISO image of Mrk 841 (first observation): the 
brightest source is Mrk 841 and the fainter one is consistent with 
PKS 1502+106; the two calibration spots are visible in the top 
and bottom left corners. 



As for Mrk 509, Mrk 841 has been extensively studied. George 

leTaLT 



1993) provided the first evidence of Fe Ka emission-fine 



detection and broad-band variability. In a sample of 24 type 1 
AGN observed by ASCA ((Reynolds] [1997}, Mrk 841 belongs 
to sources for which a "statistically significant soft excess" is 
observed. More recent observations have been performed by 

200 l|l and XMM-New ton ( |Petrucci| 
.onginotti eTaL|2010[ i: P07 studied 



BeppoSAX ( |Bianchi et al. 
[eTaL][2007l hereafter P07; 

in detail the Fe K complex variability, while Longinotti et al., 
using XMM-Newton/RGS data, investigated the warm absorber 
structure. Even though the goal of this paper is to investigate 
the broad-band X-ray spectrum of this object, in Section 3.4 we 
also discuss the Fe Ka narrow emission-line detection. 

The paper is organised as follows. In Section 2, we detail the 
Suzaku data reduction and present the serendipitous detection 
of the blazar PKS 1502+106. In Section 3.1, we study the light 
curves of Mrk 509 and Mrk 841; in Section 3.2, we analyse 
the soft X-ray excess in both sources and in Section 3.3 we 
discuss their hard excess; in Section 3.4, we concentrate on the 
Fe Ka narrow emission line detection in Mrk 841; Section 3.5 
is dedicated to the modelling of the broad-band (XIS and PIN) 
Suzaku summed spectra in a double Comptonisation scenario, 
discussing the detection of a warm absorber in both sources; 
in Sections 3.6 and 3.7, we fit the broad-band summed spectra 
with a blurred ionised reflection model and a smeared ionised 
absorption one. The discussion and the conclusions are given in 
Section 4. 



2. Data reduction 

Suzaku X-ray telescope carries two main instruments: the 
XIS ( |Koyama et aLl|2DD7] >, a soft (0.2-12 keV) X -ray imaging 
spectrometer, and the HXD flakahashi et al.|2007| l, a hard X-ray 
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detector that covers the higher energy domain (10-500 keV). 

The XIS instrument is composed of four CCD cameras: 
three of them (XIS0,2,3) are front-illuminated (FI), while one 
(XIS1) is back-illuminated (BI). Owing to an instrumental 
failure, XIS2 is unavailable for observations performed after 
November 9, 2006: for uniformity between observations, we did 
not consider data from this instrument. Event files from version 
2.0.6.13 of the Suzaku pipeline processing were used, with 
standard screening criteria, and data were then reduced using 
version 6 .7 of HEASo$\ Spectr a and light curves were extracted 
using XSelect, version 2.4a. The chosen extracting region is 
large enough (r>150 arcsec) to avoid calibration problems (as 
suggested in the Suzaku ABC Guid^J. The background was 
extracted from a region (circular or annular, according to the 
aim point) in the same image, respecting the same geometrical 
constraints as for the source. Response matrices and ancillary 
response files were generated for each XIS using xisrmfgen 
(version 2009-02-28) and xissimarfgen (version 2009-01-08) 
tools. The two XIS/FI spectra were added using mathpha 
(version 4. 1 .0.) with the 'POISS-3 ' error propagation method, 
which represents a mean value between the two extreme error 
values 'POISS-1' and 'POISS-2' (see mathpha user's guid^ 
and references therein). Response files were added using 
addrmf, version 1.21, and addarf, version 1.2.6. A single XIS 
lightcurve was obtained by summing up the three XIS0,1 and 3 
background-subtracted lightcurves, using XRonos, version 5.22. 

The HXD instrument is composed of two detectors: PIN 
(positive intrinsic negative diodes; 10-60 keV) and GSO/BGO 
(gadolinium silicate/bismuth germanate crystals; 50-500 keV). 
In this paper, we present data from PIN only, since no significant 
signal was detected from GSO. For HXD/PIN, non-X-ray instru- 
mental background (NXB) and response matrices provided by 
the HXD instrument team were used. An additional component 
for the CXB (cosmic X-ray background) ( |Boldt & Leiterj 1 987[ > 
was included and added to the NXB using mathpha. Spectra 
were corrected for the dead time, using hxddtcor, version 
1.50, and the background exposure time was increased by a 
factor of ten, following the Suzaku ABC Guide. For Mrk 841 
observations, PKS 1502+106 (see section 2.1) is in the PIN 
field of view, and its flux, as extrapolated from XIS observations, 
was added with mathpha to the background. 

Each spectrum was rebinned according to the following 
procedure: for XIS, we checked the true energy resolution of 
each instrument, extracting the spectrum of the calibration 
spots in the CCD camera (see Fig. [TJ, and measuring the 
energy and the FWHM (full width at half maximum) of the 
Mn Ka line; with these values, and using equation (1) given 
in Koyama et al. (2007), we get the true XIS energy resolution 



0.02 



0.01 



as a function of energy; for PIN, we used energy resolution 
values given in Takaha shi et al.| ( |2007[ l; we then rebinned the 
data to have five bins for each energy resolution element. To 
use a x 2 minimization algorithm, we imposed a minimum 
number of counts per bin: this second rebinning does affect 
the first and the last bins of the spectra only, while in the other 
channels the condition of five bins per resolution element is 
always dominant. We chose a minimum of 30 (for Mrk 841) and 



1 http ://heasarc . nasa. go v/lheasof t 

2 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

3 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/docs/memos/ogip 
95 008/ogip 95 008.pdf 
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Figure 2. PKS 1502+106 spectrum (observer frame) as seen by 
Suzaku/XIS/FI (black) and XIS1 (red). Data were fitted with a 
power-law model absorbed by galactic material (see section 2.1): 
data-to-model ratio is shown in the bottom plot. 



50 (for Mrk 509) counts per bin, leading to a similar number 
of bins for both sources. We considered XIS/FI data in the 
energy band 0.6-10 keV, and XIS1 data in the range 0.5-8.5 keV; 
XIS data comprised between 1.62 and 1.82 keV were ignored 
owing to a known calibration problem around the Si K edg^J 
PIN data below 15 keV were ignored, and the last PIN sig- 
nificant bin is at ^33 keV for Mrk 841 and *38 keV for Mrk 509. 

We also studied the summed spectra of the sources, adding 
XIS and PIN data from different observations. For Mrk 509, 
the first observation was not included in the summed spectrum, 
because it had been performed using the HXD aim point, instead 
of the XIS aim point used for the three other observations (see 
Table [T}. Source and background spectral files from different 
observations were added using the same procedure described 
above for the sum of XIS/FI spectra. 

To take into account the relative normalization between different 
instruments during the fitting, models for each instrument data 
were multiplied by a constant, as indicated by Maeda et al. 
(2008) Spectral fitting was performed with XSpec, version 
12.5.1. Parameter errors in both text and tables are given at the 
90% confidence level. 

We assume a flat cosmology with Ho = 70 km s^'Mpc -1 and 
Q A = 0.73. 



2.1. Contaminating source: PKS 1502+106 

Figure [Tj shows the Suzaku/XISO image of the first observation 
of Mrk 841. The brightest source at the HXD aim point position 
is Mrk 841, while a second bright source appears in the field of 
view (at a distance of about 7' from Mrk 841), that is coincident 
with PKS 1502+106. This source, situated at g = 1.839 (Smith 



|et al.|[l977) , was classified as a blazar by |Abdo et al.| (2010 



We extracted the spectrum of this source following the same 
procedure discussed in the previous section: Fig. [2] shows 
the XIS spectrum of PKS 1502+106. We fitted the spectrum 
with a simple power law absorbed by Galactic material (model 



4 http://heasarc.nasa.gov/docs/suzaku/analysis/sical.html 

5 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 



2008-06.pdf 
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wabs*powerlaw in Xspec, with Nh = 2.36 ■ 10 20 cm z ,|Dickey & 



,-2 



|Lockman| ( [T990| ). The best-fit photon index is F = 1.32 + 0.08 
iX 2 =X 2 l (degrees of freedom (DOF))=7 1 /86). The 



source 

flux in the 0.5-2 keV and 2-10 keV bands is 2.2+^ and 

7 c+0.6 
'■ J -0.7 

2+0.2 



10 13 ergs cm 2 s , 



respectively. The 2-10 keV lumi- 
nosity is 1.8^ f ■ 10 4b ergs s" 1 . In the past, PKS 1502+106 was 
observed to have a photon index F varying between 1 .4 and 1 .9 
and a 2-10 keV flux in the range 4.9 - 6.54 • 10" 13 ergs cm -2 s -1 
(|George et al.||1994| |Akiyama et al.||2003[ |Watanabe et al. 
2004) . An outburst was observed by Fermi and Swift |Abdo 
et al '' om n ^ "' ;tV - 17 - - - us m-i2<,,„„ „™-2 -i ^rr 



20101, with F . 3 _ 



0.3-i0fev = 2.18 • 10" ergs cm^ s" 1 , and 
T = 1.54. Our best-fit model parameters suggest that Suzaku 
detected PKS 1502+106 during quite a high state, but not 
during an outburst (its light curve does not show any significant 
variability). However, the power-law index is consistent with a 
low state of the source. This source falls within the PIN field of 
view, thus it has a contribution to the PIN data. The extrapolated 
flux in the 12-60 keV band is 2.5 ■ 10" 12 ergs crrT 2 s . This 
component was included in the PIN background. 
During the second observation, PKS 1502+106 falls off the 
CCD image of Mrk 841, but it remains in the PIN field of view: 
we subtracted from PIN data the same spectrum evaluated for 
the first observation. 

No bright X-ray sources are observed in the XIS field of view of 
Mrk 509. 



3. Data analysis 

3.1. Light curves 

XIS light curves of Mrk 509 and Mrk 841 are shown in Fig. [3] 
Owing to the orbit of Suzaku, the targets, periodically occulted 
by Earth, cannot be observed continuously and light curves are 
composed by several blocks of a few kiloseconds each. As no 
significant variability appears on this timescale, we decided to 
rebin the data to one bin per observation window. Light curves 
were extracted for two energy bands, a soft (0.5-3 keV) one 
and a hard (3-10 keV) one, and a count-rate hardness ratio 
(h - hard/ soft) was evaluated. The PIN light curves are always 
compatible with a constant flux. 

During the first observation of Mrk 509, the source was in 
the HXD aim point, which provides a smaller (« 10%) effec- 
tive area and a lower number of counts per second than for the 
XIS aim point (which was used for the three other observations). 
For this reason, the lightcurve of observation (obs.) 1 had to be 
corrected for this systematic error when compared to the other 
observations. Taking into account this effect, Mrk 509 displays 
a variability of as 25% between the four observations in the soft 
energy band, but within a single observation the flux variation is 
only 1 1-14%, reaching 20% only during the last observation. In 
the hard energy domain, the variability is « 16% between obser- 
vations and « 15% (but 20% during the fourth observation) on a 
10 5 s timescale. 

The value of the hardness ratio during obs. 1, 2, and 4 is 
compatible with a constant (the mean values are h = 0.27+q°,|, 
0.26 +°-°J and 0.27+jJjJ[, respectively). During the third ob- 
servation h has instead 



light curves is shown (see Fig. 3 in their paper). They found a 
correlation between the fluxes in the soft and the hard energy 
band, consistent with the nearly constancy of the hardness 
ratio, though the lower value of h during the third observation 
increases the dispersion in the correlation plot. 

The two Mrk 841 observations have an almost equal mean 
flux value, but the source shows variability on timescales of 10 5 s 
(26% and 46% in the soft X-ray domain, and 32% and 42% in 
the hard X-ray domain, for the two observations respectively). 
For this source, too, the hardness ratio is nearly constant, and 
we measure a mean value of h = 0.29^[J q 2 . 

Given the low variability observed in Mrk 509 and Mrk 841 
light curves, and to improve the signal-to-noise ratio, we study 
in the following the summed spectra as well as individual 
observations for phenomenological models. 



3.2. Soft excess 

We start the analysis considering the XIS 3-10 keV energy band. 
We model the data using a power law plus two narrow Gaussian 
functions (model zgauss in XSpec, assuming <x = 0) for the Fe 
Ka and KJ3 narrow emission lines (assuming = 7.06 keV 
and imposing its normalization factor Ckp equal to 0.15 Cxa 
( |Palmeri et al.|2003) ). An excess of emission is clearly observed 
around 6.4 keV (see Fig. [4]i and adding a broad emission line 
significantly improves the fit (F-test probability equal to 4 ■ 10" 21 
and 7 • 10" 7 for the summed Mrk 509 and Mrk 841 data set, re- 
spectively). Table[2]shows the best-fit results. The detection of a 
Fe broad emission line in Suzaku observations is consistent with 
the results showed in the FERO (Finding Extreme Relativistic 
Objects) survey ( de La Calle Perez et al.|2010 i, which studied, 
using XMM-Newton data, the relativistic Fe Ka emission lines 
in Seyfert 1 galaxies, including Mrk 509. 

Ponti etal.| ( |2009| l studied in detail the Fe K complex in Mrk 509, 
using the same Suzaku observations presented here: in particu- 
lar, they showed that a broad Fe line (with parameters consis- 
tent with our values given in Table 2) represents the best-fitting 
model, compared to a possible blend of narrow Fe lines. 
For Mrk 841, we performed the same analysis as Ponti et al. 



lower value (0.24+°°}). 



Miyazawa 



|et al.|p009] l analysed the temporal and spectral variability of a 
sample of 36 AGN observed by Suzaku: the four observations 
of Mrk 509 we present here, were studied by these authors, 
and a correlation plot between the 0.5-2 keV and the 2-10 keV ting these observations in an intermediate state). XMM-Newton 



(2009), testing the possibility that a blending of H-like and He 
like lines of Fe might be responsible for the broad line observed 
: however, this model (two narrow Gaussian lines, with energies 
fixed to 6.67 and 6.97 keV, plus the 6.4 and 7.06 Fe lines, as 
described above) does not significantly improve the goodness of 
the fit (F-test probability equal to 0.1 for the summed data-set). 
For Mrk 509, we observe variability in the 3-10 keV flux (com- 
prised between 3.39 and 3.78 ■ 10" 11 ergs cm" 2 s" 1 ), while the 
slope is consistent, within the errors, with a constant (the two 
extreme values are 1.79^SqJ and 1.86+2?^ for the second and 
the third observation, respectively). The best-fit results of the 
two Mrk 841 observations are quite similar (F = 1.75-1.77 and 
F - 1.10-1. 11-10" 11 ergs cm" 2 s" 1 ), as discussed in the previous 
section. 

It is useful to compare these results to past observa- 
tions. Using BeppoSAX, |De Rosa et al.| ((2004) detected 
Mrk 509 in two quite different states, at F2-10 keV = 2.7 and 
5.7 ■ 10" 11 ergs cm" 2 s" 1 with F = 1.80-1.59 (for compari- 
son the 2-10 keV flux of Mrk 509 during the four Suzaku ob- 



servations varies between 4.3 and 4.8 • 10 11 ergs cm 2 s set- 
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Figure 3. Mrk 509 (left) and Mrk 841 (right) light curves. Different observations are shown in different colours (one-black, two-red, 
three-green, and four-blue). From top to bottom: 0.5-3 keV light curve (in units of counts per second); 3-10 keV light curve and 
hardness ratio (3-10/0.5-3). Gap between observations is arbitrary. 
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Figure 4. Xft> data between 3 and 10 keV of the sum of Mrk 509 (left) and Mrk 841 (right) observations : data are first fitted with a 
power law (ignoring the 6-7 keV energy band, upper panel and first ratio plot); then two narrow Gaussian lines (for the Fe Ka and 
K/3 emission lines, see section 3.2) are added to the fit (second ratio plot); a broad (<x > 0) Fe emission line is finally considered 
(third and last ratio plot). Only FI data are shown, as they have the highest resolution around the Fe Ka emission line, but the fit 
was performed using both FI and XISl data. The energy on the abscissa is in the observer frame. The data were rebinned for clarity 
purposes. 



( Smit h et al.|2007[ > observed the source in a even lower state at 
F 0.5-10 keV = 2.6 • 10 -11 ergs cirT 2 s" 1 . 

Mrk 841 was observed by XMM-Newton (see P07) at 
Fi- W ke V = 0.8-1.3 ■ 10" 11 ergs cirT 2 s" 1 with F = 1.30- 
1.95. Compared to these values, the source was in an 



intermediate state during Suzaku observations; in the past, 



Mrk 841 was seen in a quite low state by Pounds et al. 



([1994b (simultaneous Ginga and ROSAT observations) with 
^2-10 keV = 7.9 ■ 10~ 12 ergs cm" 2 s -1 . 
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Table 2. Best-fit parameter values for the model powerlaw + zgauss + zgauss + zgauss, considering XIS data between 3 and 10 
keV. The three Gaussian functions model the Fe Ka and K/3 narrow emission lines (with cr = 0) and the broad Fe K emission line. 
The Fe K/3 emission line was fixed in energy (Exp = 7.06 keV) and in normalization (imposed equal to 0.15 times the Ka line 
normalization). F is the power-law photon index; the power-law normalization Cp is the number of photons keV cnT 2 s _1 at 1 
keV; the line normalization parameters Ck» and Cbr are the integrated number of photons cirT 2 s _1 in the line; <x ir is the line width 
expressed in keV; EqW is the line equivalent width expressed in eV; the flux F, evaluated between 3 and 10 keV, is given in units 
of ergs cirT 2 s _1 . For the summed spectra, the results of a XIS (above 3 keV) + PIN fit are reported as well (last line): the R pexrav 
parameter represents the pexrav reflection scaling factor (see Section 3.3). 
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Figure 5. Mrk 509 (left) and Mrk 841 (right) soft excess: data from 3 to 10 keV are fitted with a power law plus three Gaussian 
emission lines (see Table [2]), and the model is then extrapolated to lower energies. Galactic absorption is taken into account. Only 
XIS1 data are shown, as they extend to 0.5 keV, but the analysis was performed using both FI and XISl data; the colour code for 
the different observations is as in Fig. [3] The energy on the abscissa is in the observer frame. The data were rebinned for clarity 
purposes. 



Extrapolating the power-law fit between 3 and 10 keV to 
low energies, an excess of photons is observed: in Fig. [5] we 
report XISl spectra for Mrk 509 and Mrk 841, and for each 
observation; the bottom panel shows the data-to-model ratio. 
The Galactic absorpt i on has been taken into account following 
Dickey & Lockman| ( |1990| , using N H = 4.11 ■ 10 20 cirT 2 for 
Mrk 509 and N H = 2.34 • 10 20 cnr 2 for Mrk 841. The soft excess 
shape is remarkably similar for the two Seyfert galaxies, and for 



Mrk 509 spectral variability is observed at low energies. 
As a first qualitative study, we model the 0.5-10 keV spectrum 
with a broken power law (plus three Gaussian emission lines, 
with parameters fixed as in Table [2] and Galactic absorption, 
fixed at values given above) to evaluate the soft excess slope 
(r s ) and energy break Ef,. As shown in Table [3] this model does 
not provide a good fit to our data : the soft spectrum between 
0.5 and 3 keV clearly displays some curvature and cannot be 
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described by a power law. 



3.3. Hard excess 

Extrapolating the 3-10 keV best-fit power law to high energies, 
using PIN data, a hard excess can also be seen. This feature 
is well-explained by introducing a distant reflection by cold 
material, which accounts at the same time for the Fe Ka 
narrow emission line and the high energy bump. To quantify the 
amount of reflected emission required, we fitted the summed 
PIN and XIS (above 3 keV) data with a power law (plus three 
Gaussian lines, as described in the previous section) plus pexrav 
(Magdziarz & Zdziarski 1995 ), a model that describes reflection 



from neutral material. The pexrav free parameters are the 
reflection scaling factor R peX rav (which is equal to 1 for a source 
isotropically illuminating the disc); the illuminating power-law 
high-energy cut-off E c (fixed at 100 keV); the abundances (fixed 
to equal solar values) and the disc inclination angle 8 (fixed at 
30 degrees). For the summed spectra, we obtain R pexrav — 0.4^2 

and 1.3+Q3 for Mrk ^ 09 and Mrk 841 > respectively. Adding 
the neutral reflection component slightly modifies the best-fit 
parameters of the broad Fe emission line, in particular its equiv- 
alent width. Best-fit parameters are given in Table 2. Pexrav 
has the advantage that gives directly the information about the 
reflected fraction R 



pexrav: 



but does not evaluate the narrow Fe 
Ka emission line produced. We then decided to use reflext, a 
modified (and not public) version of reflion model (an older 
version of reflionx (Ros s et al.|1999||Ross & Fabian|2005] l; see 
section 3.6). Unlike reflion, in which the ionisation parameter 
£ is comprised between 10 and 10 4 erg cm s , in reflext the £ 
parameter space was extended as low as 1 erg cms -1 . Using 
reflext, the evaluation of the reflected fraction is less accessible, 
but the Fe Ka narrow emission line and the high energy bump 
are intrinsically tied. 

In the following, we add reflext with parameters FelFe Q = 1, 
r 'reflext = P, and £ = 1 erg cm s . Its normalization parameter 
Crefiext was fixed to ensure that the flux of the model in the 
narrow Fe Ka emission line reproduced the observed line flux. 
As the intensity of the reflected emission depends on the slope 
of the illuminating continuum, we evaluated the value of C re fi ext 
that reproduces the observed line flux as a function of T; we 
fitted the result with a polynomial function, and we fixed the 
reflext normalization as a function of T. 



Table 3. Best-fit parameter values for the broken power-law 
model. XIS/FI data between 0.6 and 10 keV and XIS1 data be- 
tween 0.5 and 8.5 keV were considered. The power-law normal- 
ization Cu is the number of photons keV -1 cm -2 s -1 at 1 keV. 
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Figure 6. Contour plots (in black for obs. 1, and in red for obs. 2) 
of line flux versus line energy (source frame) for the narrow 
emission Gaussian line in Mrk 841. The three confidence levels 
are for 68, 90, and 99%, respectively. In grey, we superimpose 
the contour plots of 2001-2005 XMM/Newton observations (see 
P07, Fig. 9). 



iron Ka for both observations. 



3.4. Fe Ka narrow emission line in Mrk 841 

We present here the analysis of the Fe Ka emission line in 
Mrk 841. Ponti et al. (2009 1 already presented a thorough study 
of the iron K complex in Mrk 509, using the same Suzaku 
observations we analysed. 

A narrow (cr = 0) Fe Ka emission line was observed in previous 
Mrk 841 observations. In particular, P07 claimed that this 
narrow component is roughly constant in time, in agreement 
with a remote reflection scenario. We tested this hypothesis 
fitting the 3-10 keV XIS data with a power law plus a narrow 
Gaussian emission line, with cr fixed to zero, as done by P07. In 
Fig.|6| we report the contour plots of line flux versus line energy, 
for the two Suzaku observations, and we compare them to the 
results of P07 (XMM/Newton, observations from 2001 to 2005). 
The line flux agrees with the findings of P07, and no variability 
is observed. Line energy, in the source frame, is consistent with 



In the following, we concentrate on the broad-band spectral 
fitting of Mrk 509 and Mrk 841 data with three different physical 
models : a double Comptonisation, a blurred reflection, and a 
smeared absorption. As only weak flux variability is observed 
between individual Mrk 509 observations, and as no significant 
variability is observed between individual Mrk 841 observations, 
we perform and discuss the broad-band spectral analysis for the 
summed spectra only (unless explicitly mentioned). 

3.5. Double Comptonisation model 



In the standard disc-corona scenario ( Haar dt~& Maraschi| 
|1993 1, the power-law high-energy emission is explained as 
a Comptonisation of thermal photons emitted by the black 
hole accretion disc. Electrons producing this inverse Compton 
effect have a temperature kF e * 100 keV and a low optical 
depth, t sis 1, One possible way of easily explaining the soft 
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excess is that there is a second warm Comptonisation region, 
of lower temperature and higher optical depth, that modifies 
the black-body emission from the accretion disc, scattering UV 
photons up to soft X-ray energies. A physical interpretation of 
this "second corona" might be a warm skin on the disc surface 
( Janiuk et al.|200l] >, or a transition region between an outer disc 
and an optically thin inner flow (Magdzia rz et al.| [l998). We 
used NTHComp (Zycki et al. 1999 1, which is one of the XSpec 
models, to simulate Comptonised spectra: this model solves the 
Kompaneets equation, and the output spectrum is parametrised 
by the electron temperature kT e and the spectral index T; seed 
photons come from an accretion disc spectrum, parametrised 
by kTdisc, the temperature at the inner radius. The electron 
temperature of the corona producing the high energy emission 
(subscript h), IcT/,, can be constrained only via the high energy 
cut-off: since we do not observe this cut-off in our data, we fixed 
it at 100 keV. For the soft Comptonisation region (subscript s), 
we cannot fit at the same time F s and kT s , and, following GD04, 
we set F s = 2. 

The exact value of kT^ sc does not affect the fit as long as it 
is lower than our energy domain, thus we chose kTdi SC = 45 
eV. We added reflext to our model to describe the reflection 
from distant material, with a normalization parameter fixed as 
described above, plus two Gaussian functions for the Fe K/3 and 
the Fe Kj, roa d emission lines, as in Table [2] (for the XIS + PIN 
fit). Our best-fit model results are given in Table [4] 

The double Comptonisation model was studied in detail by 
GD04, who fitted a sample of 26 AGN showing that the cold 
Comptonisation region has a remarkably constant temperature 
<kT s >- 0.12 keV with cr = 0.02 keV, despite large variations in 
black hole mass and AGN luminosity. They also used the model 
NTHComp. The values we found for Mrk 509 and Mrk 841 
during Suzaku observations agree with their statistics. We recall 
that the exact value of kT s depends on the assumption made for 
r s : the only physical parameter that can be constrained (when 
the high energy cut-off is not observed) is the Comptonisation 
parameter y, defined as 4kT/(mc 2 ) ■ Mox(t,t 2 ), where r is the 
plasma optical depth. The T parameter represents the asymptotic 
power-law function that describes the Comptonised spectrum, 
and can be expressed as (2.25 + 4/y) 5 - 0.5, which is equal to 
2 for y — 1 . The most important result obtained by GD04 is not 
the exact value of kT s , but its constancy in AGN samples. 
Following GD04, we evaluated the soft-excess strength, Rcd, 
defined as the ratio of the unabsorbed 0.3-2 keV flux in the 
cool to that in the hot Comptonised components: it is equal 
to 0.23+°$ for Mrk 509 and 0.28+°;°^ for Mrk 841, placing 
our sources in the lowest part of the GD04 sample in terms of 
soft-excess strength. 



3.5.1. Additional contribution of the warm absorber 

In terms of the y 2 values, the double Comptonisation model 
does not provide a good fit to the data. For both sources, the 
presence of a warm absorber was confirmed, and its features 
could explain the residuals of the fit. The most complete studies 
of the warm absorber of these two objects were done using 
XMM/RGS data by |Smith et aTJ(|2W7) and|Detmers et al.| ( [20T0l > 



Longi notti etaLf ( |2010| l for Mrk 841 



for Mrk 509, and by 

Smith et al. detected a three-phase warm absorber in Mrk 509, 
consisting of a low ionisation component (log £ * 0.89 and 
column density Nh ~ 7.9 • 10 20 cm -2 ), an intermediate ionisation 
phase (log £ m 2.14, Nh ~ 7.5 • 10 20 cirT 2 ), and a high ionisation 



one (log £ ss 3.26, Nh ~ 55 ■ 10 20 cirT 2 ). Detmers et al. 
confirmed the presence of three different components in the 
warm absorber: they found a low ionisation phase (log £ * 0.6, 
Nh ~ 1.0 • 10 20 cnT 2 ), an intermediate ionisation phase 
(log £ * 1.95, N H ~ 10.5 • 10 20 cm~ 2 ), and a high ionisation 
phase (log £ ss 3.2, Nh ~ 80.0 ■ 10 20 cirT 2 ). 
A multi-phase warm absorber was also found in Mrk 841, by 
Longinotti et al., with an intermediate ionisation component 
(log £ * 1.5-2.2, Nh ~ 12-39 ■ 10 20 cirT 2 ) and a higher ionisation 
phase (log f = 2.8-3.3 and Nh ~ 76-300 ■ 10 20 cnT 2 ). 
We tested the presence of the warm absorber using the model 
zxipcf ( |Reeves et al.||200"8) , which is another of the XSpec 
models, which has as free parameters the ionisation parameter 
log the column density A^, the covering factor /, and the 
redshift z. We started our study by multiplying the double 
Comptonisation model by one zxipcf component, with / = 1 
and z = Zobject- For Mrk 509, we found that log % WA = 2.19+° ^ 
and Nh,wa = 17+J • 10 20 cirT 2 , values that are consistent with 
the intermediate ionisation phase found by Smith et al. and 
Detmers et al.. For Mrk 841, we found that log g WA = 2.3+°,J 
and Nh.wa = 17+ 4 ■ 10 20 cirT 2 , which is consistent with the 
intermediate ionisation phase in Longinotti et al.. Adding zxipcf 
always improves the goodness of the fit, as is confirmed by 
a F-test: the null hypothesis probabilities were found to be 
7 ■ 10~ 44 (for Mrk 509) and 1 • 10~ 13 (for Mrk 841). 
We then tried to add a second absorber (subscript 2). For 
Mrk 509, we found that a two-phase warm absorber is statis- 
tically preferred (the associated F-test probability is 5 ■ 10~ 5 ) 
and we measured log %wa,2 = 2.7+° ) j, N h ,wa,2 = 7+? " 10 20 cm~ 2 , 
log £ WAA = 1.7+g;}, and N h ,wa,\ = 8! 2 ■ 10 20 cnT 2 . For the high 
ionisation phase, XMM/RGS data also showed a warm-absorber 
flow velocity blueshifted relative to the galaxy: we then fixed 
Zwa,2 to the value found by Detmers et al. (zwa,2 = 0.0334, 
corresponding to v = -290 km s^ 1 ). For M"k 841, adding a 
second absorber does not improve the goodness of the fit and 
the absorber parameters are poorly constrained: we then fixed 
log %wa,2 = 3.3, which is equal to the value measured by 
Longinotti et al. during the 2001 XMM observation of Mrk 841. 
Our best-fit results are given in Table |4] and plotted in Figs. [7] 
and|8](left panel). We considered whether there is a three-phase 
warm absorber in Mrk 509, but the goodness of the fit worsens, 
even fixing the warm absorbers parameters at the values mea- 
sured by XMM/RGS. 

The Zxipcf model can describe partial covering absorbers, 
with /<1. We tested the effect of this parameter on a one -phase 
warm absorber. For Mrk 509, leaving / free to vary, it decreases 
to 0.3, while Nh.wa increases, reaching 3 ■ 10 22 cirT 2 ; the new fit 
is statistically better than the one with / = 1 (the null hypothesis 
probability is equal to 1 ■ 10~ 7 ). For Mrk 841, the best-fit result 
is always consistent with / = 1 (we measure / > 0.7). Despite 
the better fit obtained with f<\ for Mrk 509, we decided to 
use and study the / = 1 case in order to compare our results 
to those obtained with the XMM/RGS data (and for consistency 
between the analysis of our two sources). In any case, there is 
a degeneracy between / and Nh, and leaving / free affects the 
values of £ and T s and, consequently, the physical interpretation 
of our results. In the next sections, we always use / = 1, with a 
reminder that, for Mrk 509, allowing / free to vary always leads 
to a better fit. 

In both sources, XMM/RGS data revealed an OVII emission 
line at E ss 0.56 keV. We tested for the presence of this emission 
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by adding a narrow (<x = 0) Gaussian line to the model (the fit 
was performed using data below 10 keV, without reflext, which 
produces an OVII emission line as well (see below)). 
For Mrk 509, we obtained a significantly better fit (F-test 
probability = 4 ■ 10" 9 ) for the fourth observation, while for 
the other observations and the summed spectra the goodness 
of the fit did not improve. However, at 0.56 keV we do not 
have the XIS/FI detectors (which start at 0.6 keV), and we are 
at the energy edge of XIS1: the energy resolution is quite low 
(the emission line is defined by one bin only) and we cannot 
confirm whether this statistically significant improvement to the 
fit corresponds to a true emission-line detection. 
For Mrk 841, adding a Gaussian OVII emission line does 
not improve the goodness of the fit in either any individual 
observation or the summed spectrum. 

The origin of this emission line is thought to be the warm 
gas itself but we do also expect OVII emission from colder 



reflection. In particular, Ebrero et al. ( 2010) > showed, for the 
Seyfert 1 Mrk 279, that the reflext code (for log £ < 1.43) 
produces an OVII line emission flux that exceeds the observed 
OVII line flux by a factor of 10 or more. For Mrk 509 and 
Mrk 841, this excess does not exist: the line flux produced 
by reflext (with log £ = 0) is equal to 2.1 • 10~ 4 cirT 2 s for 
Mrk 509 (compared to 2.7 • 10~ 4 cirT 2 s _1 found by Detmers 
et al.), and to 7.2 • 10~ 5 cm" 2 s _1 for Mrk 841 (compared to 
(6.4 ± 5.3) ■ 10~ 5 cirT 2 s _I found by Longinotti et al. in the 
January 2005 XMM/RGS observation). 

GD04 performed their fits without testing for the presence 
of a warm absorber in their data. Looking at Table|4] we can see 
that, after including zxipcf, kT s increases from 0.13-0.14 keV to 
0.19-0.21 keV, becoming only marginally consistent with GD04 
results. 



3.6. Blurred reflection model 

The constancy and universality of the soft excess temperature 
is difficult to explain in a double Comptonisation scenario, for 
objects with different black hole masses (GD04). It is instead 
expected if the soft excess origin is linked to atomic processes. 
If the reflection from the accretion disc is ionised, and if the 
emission lines are broadened by a relativistic velocity near 
the black hole, this component might also explain the soft 



excess. We used the reflionx model ( |Ross et aL]|1999| Ross &| 
Fabia n|2005| l, convolved with the kdblur model ( |Laor|1991| ), to 
describe the disc ionised reflection. For the relativistic blurring, 
we assumed that R/„ = 1.235 Rc (the last stable orbit for 
a near extremal Kerr black hole with parameter a = 0.998, 
( |Thome|[T9"74j l), R out = 100 R G , = 30 degrees, and disc 



emissivity index f$ — 3 (|Laor 1991 ); for reflionx, we fixed 



the iron abundance to equal the solar value, leaving £ and the 
normalization free. Following what we did in the previous 
section, we added a neutral reflection from distant material 
(reflext, with normalization parameter fixed as a function of T, as 
described in Section 3.3), a narrow Gaussian line (for the Fe K/3 
emission line), and a two-phase warm absorber totally covering 
the source (zxipcf), with an ionisation parameter and column 
density left free to vary within the error bars of the values found 
using a double comptonisation continuum (see Section 3.5.1). 
A broad emission feature around 6.4 keV is expected from the 
blurred reflection on the accretion disc, so no broad emission 
line was included in the model. Best-fit model results are given 
in Table |5] 



The R parameter, which gives information about the reflected 
flux fraction, as defined by C06, is the ratio of the unabsorbed 
reflected component (reflionx) flux to the total (excluding 
reflext) unabsorbed flux in the 0.3-12 keV energy band. For a 
flat disc illuminated by an isotropic source, we evaluated a value 
of R a 0.5 for a high ionisation state (log £ > 3.5), while fol- 
lower ionisations the reflection fraction decreases, reaching 0.4 
for log £ = 3, 0.2 for log £ = 2, and 0. 1 for log £ = 1 .5. 
This model is widely used in the literature: C06 analysed a 
sample of 34 type 1 AGN, fitting their XMM-Newton obser- 
vations with the reflion model (an older version of reflionx). 
For the R parameter, they found that one third of the sources 
have no power-law component (R = 1), while in the rest of 
the sample R is almost constantly distributed between 0.25 
and 0.8. Our values, which range between 0.15 and 0.18, are 
quite low compared to this result. In particular, Mrk 841 is 
included in their sample, and they found that R = 0.67 + 0.17. 
This value differs from what was found by P07, in a dedicated 
paper on XMM-Newton observations, including the one studied 
by C06, for which they found that R = 0.3 (in their paper, 
P07 evaluated that R' = 0.3, where R' is the reflected fraction 
evaluated in the energy band 0.1-1000 keV, when including 
warm absorption; we reevaluated R using their best-fit model 
parameters). In our model, we froze all kdblur parameters. We 
tried to allow /3 to vary, but its best-fit value is always consistent 
with 3. C06 performed their fits leaving all, apart from R out , of 
the kdblur parameters free to vary: with such a large number of 
parameters, best-fit model values become less well-constrained, 
with several parameters pegged at either lower or higher limits 
(as can be seen in C06, Table 3) and a physical interpretation 
of these results can be misleading. The results of P07 for their 
XMM-Newton observations of Mrk 841 were 1.59 < F < 2.45, 
2 < log £ < 2.5, and /3 ~ 4 (but in one observation they found 
that /3 > 8.6), with a reflection fraction R' varying between 0.30 
and 0.43. 



The best-fit model result for Mrk 509 is not good (x 2 = 1.25) 
and some residuals are clearly observed between 6 and 7 keV 
(see Fig. [9 1. We first tried to vary the inner radius and the incli- 
nation angle of the accretion disc, significantly improving the fit 
(F-test probability = 1 ■ 10~ 4 ) and measuring R, n < 1.7 Rc and 
8 - 38^ 2 degrees. However, this cannot take into account all the 
residuals observed near the Fe Ka line. An excess of emission in 
the Fe Ka complex might be caused by an iron abundance higher 
than the solar value: we tried to vary this parameter but the best- 
fit model results are always consistent with Fe/Fe = 1. Adding 
a broad Gaussian line, the fit is significantly improved (the null 



hypothesis probability is 10 ). We measured Eoj, roat / = 6.9 



+0.1 
-0.1 

keV, C broad = 10+] • 10~ 5 photons cirT 2 s" 1 , and cr hroad = 0.8^ 
keV, with an equivalent width of 255+gg eV (see Table |5] and 
Figj7] middle panel). 

For Mrk 841, a good fit to the data was achieved using the first 
model (all kdblur parameters were kept fixed at the values de- 
scribed above, and there was assumed to be no broad emission 
line), but the fit improved when 6 was allowed free to vary, and 
we found that 6 — 41^ degrees (with a F-test probability equal 
to 4 • 10~ 4 ). Adding a broad emission line at the Fe K complex 
energy does not improve the goodness of the fit (F-test proba- 
bility = 0.1 compared to the fit where 9 is allowed to vary) (see 
Table [5] and Fig j5] middle panel). 
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Figure 9. Ratio of the Mrk 509 data to the blurred reflection 
model either with (bottom) or without (top) a broad emission Ka 
line (see section 3.6). XIS/FI data are in black, while HXD/PIN 
data are in green. For clarity, the data were rebinned, and XIS I 
data are not shown in this plot. 



3. 7. Smeared absorption model 

An alternative model to explain the soft excess is a relativisti- 
cally smeared absorption, where high velocities can broaden 
all atomic features, producing a curved continuum spectrum. 
This absorption might correspond to an outflowing disc wind, 
as shown by |Murray & Chiang] ( |1997| >. We used swind (GD04), 
one of the XSpec models, to describe this component: a XSTAR 
(Bautista & Kallman 2001) grid (calculated assuming T = 2.4 
and a microturbulent velocity of 100 km/s) is convolved with a 
Gaussian function of width <x = v/c. Free parameters of swind 
model are the column density Nh, the ionisation parameter 
£ and the velocity dispersion cr. As for the ionised reflection 
model, we added a neutral reflection (reflext), a narrow Gaussian 
line for the Fe KB emission line and a two-phase warm absorber 
(zxipcf model). As no broad Fe emission line is included in the 
model, we also added a broad emission line with parameters 
fixed as in Table|2](for the XIS + PIN fit). Best-fit model results 
are shown in Tableland Figs. [7] and [8] right panel. 
For both sources, the values of cr are always pegged at the 
highest authorized value (v = 0.5 c), and we can only give lower 
limits. It is interesting to cite the results of a fitting without 
zxipcf: cr decreases to * 0.3 (it is equal to 0.29+°jjJ, and 0.28+°^ 
for Mrk 509 and Mrk 841, respectively). We can compare 
our results for Mrk 841 with P07, who tested the same model 
for XMM-Newton observations: their best-fit values for swind 
model are 2.4 < log £ < 3.5, 0.9 < N H < 4.4 • 10 23 cirT 2 , and 
cr consistent with 0.3 ( but > 0.4 in two cases). The values of 
Nh and £ that we found agree with these results. This model 
was studied by |Middleton et al.| ( |2007| ) (hereafter M07), who 
analysed the same AGN sample used by C06. They found that, 
in half of their sample, cr is pegged at 0.5, and that, in general, 
it has a typical value > 0.3. Their £ and Nh best-fit values do 
not have large dispersion: log £ varies between 2.5 and 3.7, with 
lower values (log £ < 3) for their narrow-line Seyfert 1 (NLS1) 
galaxy sample; the values of Nh are between 5 -10 22 and 5- 10 24 , 
with, once again, values lower than 1.5 • 10 23 restricted to NLS1 
galaxies. The intrinsic power-law spectral index varies between 
2.0 and 3.4, with a clear difference between the softer NLSls 
and the harder broad-line AGN. The values that we found agree 
with these global results, but F is only marginally compatible 
with M07's results. The evaluation of F is important since it 



affects the soft excess strength and shape, and therefore the 
swind parameters. To examine whether PIN data have an effect 
on the evaluation of F, we tried to fit XIS data alone, but the 
best-fit values of F are compatible with the broad-band fit. 



4. Discussion and conclusions 

Both Mrk 509 and Mrk 841 are bright Seyfert 1 galaxies 
that are well-known for their soft X-ray excess, their iron Ka 
emission line, and their spectral variability. We have presented 
our analysis of new Suzaku observations: four of Mrk 509 
(from April, 2006 to November, 2006), and two of Mrk 841 
(January and July, 2007). The soft excess has been found to be 
clearly present in both sources, with a flux as 2 times higher 
than the 3-10 keV power-law extrapolation at 0.5 keV. |Pounds| 
& Reeves (2002) analysed a sample of six Seyfert 1 galaxies, 
showing that a variety of soft excess shapes exists, correlated 
with AGN luminosity: our sources during Suzaku observations 
can be described as 'gradual soft excess' galaxies, according to 
Pounds and Reeves terminology. We can compare our results 
to pas t observations (pounds et al. 200 1| Pounds & Reeves 



2002; |De Rosa et aL||2004[ P07): a soft excess two times 
higher than the power-law extrapolation at 0.5 keV is common 
for both sources, even though Mrk 841 during XMM-Newton 
observations showed a stronger soft excess, 2.5-3 times higher 
than the power law, and the soft excess of Mrk 509 during obs.l 
is quite low, similar to what was observed during the highest 
state in BeppoSAX observations ( |De Rosa et al.|2004] ). 

Comparing our observations to statistical studies of the 
soft excess (GD04, C06, M07), we have confirmed that both 
sources were detected during a weak soft-excess state, setting 
them always in the lower part of the sample: for the double 
Comptonisation model, the "soft excess strength" (Rcd, as 
defined by GD04) of Mrk 509 and Mrk 841 during our obser- 
vations is between 0.23 and 0.28, while the GD04 sample has a 
typical value of 0.4-0.6; in the reflection model, we found a R 
parameter between 0.15 and 0.18, which is quite low relative to 
the reflected flux fractions of C06. 

Both sources exhibit a hard excess at energies greater than 
10 keV. This excess is weak compared to what was observed by 
Suzaku in three other Seyfert 1 galaxies (Wal ton et al.||2010] l, 
though they are rather extreme exemples of hard excesses. 

A narrow Fe Ka emission line is detected in both sources, 
consistent with a reflection by remote and cold material (i.e. 
the dust torus or the furthest regions of the accretion disc). For 
Mrk 841, we compared our results to previous XMM-Newton 
observations, showing that the line flux is constant in time. For 
an analysis of the Fe K complex emission of Mrk 509, we refer 
to |Ponti et al.| ( |2009] ). A broad Fe K emission line is observed as 
well, with equivalent width ss 150 eV in both sources. 

We detected a two-phase warm absorber in Mrk 509 data, 
with log £wa,\ ~ 1-7 and log %wa,2 ~ 2.7, values that are 
lower than those found by Smith et al. and Detmers et al.. In 
the Suzaku data, we did not detect the low ionisation phase 
(log £, = 0.6 - 0.9) of the absorber observed by XMM/RGS. 
For Mrk 841, we found a low-ionisation warm absorber with 
log £wa,i ~ 2-2, consistent with the lower ionisation phase in the 
two-phase absorber described by |Longinotti et aT (2010 1, and 
we fixed the high ionisation parameter to be log £wa,2 = 3.3. For 



M. Cerruti et al.: Suzaku broad-band observations of the Seyfert 1 galaxies Mrk 509 and Mrk 841 



11 



Mrk 509, considering a warm absorber covering factor lower 
than 1, always improves the fit. 

We started the analysis of the broad-band spectra of the 
two sources by considering the double Comptonisation model, 
which correctly fits the data (x 2 w 1.09 for Mrk 509, and 
X 2 ~ 1.00 for Mrk 841), but the constancy of the temperature of 
the cold Comptonisation region in statistical studies of Seyfert 1 
galaxies is not explained by current physical models of accretion 
disc atmospheres, and thus disfavours this model. 
We then fitted the broad-band spectrum with the two compet- 
itive models developed to explain the soft excess with atomic 
processes: a blurred reflection and a smeared ionised absorption. 
For M"k 509, we need to add a broad emission line to correctly 
fit the data, and we obtain y 2 ~ 1.04-1.02 for the blurred 
reflection and the smeared absorption model, respectively. 
The physical origin of this emission is unclear: for the reflection 
model, a broadened emission line is already included in the 
code we used, and it is difficult to explain why we do need some 
more reflection at E « 6.4 keV. Interestingly, this result suggests 
that the reflection model correctly describes the soft excess 
shape and strength, but underestimates the broad emission at Fe 
K complex energies. 

For the smeared absorption model, we did not include any 
reflection from the accretion disc, and the fact that we need 
to add a broad emission line is less critical: the swind code 
assumes a simple power-law emission as input for the smeared 
absorption, and the presence of residuals in the Fe K complex 
region means that the continuum spectrum absorbed by the 
wind does include, indeed, some reflected component. 
For Mrk 841, no extra emission in the reflection model is 
required, while a broad Gaussian emission line is needed in the 
absorption model. Each model gives a good fit to the data, and 
none of them can be invalidated (x 2 as 1.04 for the two models). 

Suzaku data, characterised by a good energy resolution 
in a large energy domain (0.5-60 keV), seemed to be able to 
constrain the intrinsic power-law photon index at high energies, 
and then to favor one model over another. The best-fit models 
show that, while in the soft X-ray domain the models are 
degenerate, in the hard energy band they do not superpose any 
more and we could, in principle, distinguish between them, but 
the PIN signal-to-noise ratio is still too low to constrain the fit. 
For the sum of Mrk 509 observations, we can estimate from the 
best-fit model in the hard X-ray domain a flux F n-60fev - 4.47, 
4.81, and 4.27 • 10" 11 ergs cirT 2 s" 1 and F 14 _ 150teV = 7.64, 8.07, 
and 7.66 • 10~ n ergs cm~ 2 s for the double Comptonisation, 
the blurred reflection, and the smeared absorption model, 
respectively; for the sum of Mrk 841 observations we obtain 
Fn-60kev = 1-94, 2.15, and 1.67 ■ 10~ n ergs cirT 2 s" 1 and 
Fu-wkeV = 3.47, 3.72, and 3.06 ■ 10" 11 ergs cirT 2 s" 1 for the 
three models, respectively. These values can be compared with 
the INTEGRA L ( |Krivonosetal.|2010| l and Swift/BAT ( |Cusumano| 
|et al.|[20T0) i detections of our sources: INTEGRAL detected 
Mrk 509 at F 17 _ 60 k e v = 5.90- 10" 11 ergs cirT 2 s _1 and Mrk 841 at 
Fn-60kev = 3.07 ■ 10 _1 1 ergs cirT 2 s _1 ; Swift/BAT observed a flux 
in the 14-150 keV energy band equal to 8.5 • 10" ergs cirT 2 s _I 
for Mrk 509 and 2.9 ■ 10" ergs cirT 2 s" 1 for Mrk 841. The 
differences between the predicted high-energy fluxes demon- 
strate that a simultaneous broad-band observation with adequate 
signal-to-noise ratio would be able to distinguish between the 
different models. 

Statistical studies of AGN can help us to solve the problem of 



the origin of the soft excess. The discovery of the remarkably 
constant temperature of the cold Comptonisation region, despite 
the AGN luminosity and black hole mass spanning several 
decades, is hardly explained in a double Comptonisation 
scenario. Our kT s values agree with GD04 results, but we note 
that adding a warm absorber (which was not included in the 
GD04 analysis) affects the kT s best-fit values. 
Reflection and absorption models, that are tied to atomic 
processes, can easily explain the small spread in the soft-excess 
temperature. However, the smeared absorption model has 
the weakness that the velocity dispersion required to fit the 
soft excess is very large, reaching c/2, which corresponds to 
extremely relativistic terminal velocities of the wind. Numerical 
simulations ( |Schurch & Done| |2008| l show that disc winds 
have cr < 0.1c, which is very low compared to swind model 
results. An outlet could be a magnetically driven wind, partially 
covering the source. We underline that swind takes into account 
only the absorption by the wind, but if we also consider its 
emission, ( [Chevailier et al.|2006|rSchurch & Done|2006| l the fit 
requires a lower smearing velocity and a lower column density. 
A possible concern of the swind model is that the absorption 
is evaluated assuming an incident power-law spectrum with 
r = 2.4, which is inconsistent with the F best-fit values of our 
observations. 

The Fe broad-line excess we observe in both sources suggests 
that, even in an absorption scenario, we need some reflected 
emission from the inner regions of the accretion disc. A broad 
Fe K emission line can also be reproduced by disk wind models 
( |Sim et al.|2010[ ), with lower outflow velocities (0.1c). 
The blurred reflection model fits the Suzaku data well for 
Mrk 841, while for Mrk 509 we need to include an additional 
emission at Fe K energies, which is difficult to explain. The 
model does not need a high reflection fraction R, which is seen 
in several other cases, and is explained by a huge light bending 
effect ( IMiniutti & Fabian|2004| close to the black hole. 
It is of course possible that different processes are at work 
at the same time. In this scenario (which is difficult to test 
in a stationary framework owing to the high number of free 
parameters involved), the broad emission observed in Mrk 509 
could be due to standard ionised reflection from the accretion 
disc, while Comptonisation or smeared absorption should be 
responsible for the soft excess below 1 keV (see P atrick et al.| 
POTT) ). 

Alternatively, the continuum X-ray emission from the hot 
corona might deviate from a power law when measured over a 
broad energy band, showing an intrinsic curvature that has not 
been taken into account in our modeling. 

Variability studies can help us ascertain the physical origin of 
the soft excess. In the blurred reflection scenario, variability 
can be explained by changes in the height of the corona above 
the disc, but this picture implies a strong continuum-reflected 
component correlation that has not been observed. In the 
smeared absorption model, spectral variability is more easily 
explained by a variation in the absorber physical parameters, or 
its covering factor. Unluckily, neither of our sources displayed 
significant variability during Suzaku observations, and we could 
not use temporal information to more tightly constrain our 
models. 



Future X-ray observatories, and in particular Astro-H ( Takahashi 
|etal.|2010| > and NuSTAR ( |Harrison etal.|2010) , wich have higher 
sensitivities and energy resolutions in the hard X-ray domain, 
will be able to more tightly constrain the spectrum, casting light 
on the physical origin of the soft excess in Seyfert 1 galaxies. 
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Table 4. Best-fit parameter values for the double Comptonisation model, with (second row) or without (first row) a two-phase 
warm absorber. Fs and IcTh are fixed at values 2 and 100 keV, respectively. Distant reflection is modelled using reflext, whose 
normalisation parameter is fixed to reproduce the observed Fe Ka line (see Section 3.3). Two Gaussian functions are included in 
the model, for the narrow Fe K/3 and the broad Fe K emission line, with parameters fixed as in Table [2] NTHComp normalization 
parameters C are in units of keV cirT 2 s _1 at 1 keV; warm absorber column density Nr is given in units of cm -2 , and ionisation 
parameter £ in units of erg cm s . Fixed parameters are marked by (F). The last two columns give the ratio of the chi-squared value 
to the number of degrees of freedom and the associated probability. 
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Table 5. Best-fit parameter values for the blurred reflection model. Distant reflection is modelled using reflext, whose normalisation 
parameter is fixed to reproduce the observed narrow Fe Ka line (see Section 3.3). A Gaussian function is included in the model, for 
the narrow Fe K/3 emission line, with parameters fixed as in Table [2] For Mrk 509, the best-fit results including a broad Gaussian 
emission line are also shown, reflionx and power-law normalization parameters C are in units of keV~' cirT 2 s _1 at 1 keV; ionisation 
parameter £ is given in units of erg cm s -1 ; the line normalization parameter C 'broad is the integrated number of photons cirT 2 s~' in 
the line; cr is the line width expressed in keV. The R parameter is defined as the ratio of the flux of the unabsorbed reflected com- 
ponent to the total unabsorbed flux (evaluated between 0.3 and 12 keV). The two-phase warm absorber parameters are constrained 
within the range evaluated with the double Comptonisation model (Table [4]). The last two columns give the ratio of the chi-squared 
value to the number of degrees of freedom and the associated probability. 
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Table 6. Best-fit parameter values for the smeared absorption model. Distant reflection is modelled using reflext, whose normal- 
isation parameter is fixed to reproduce the observed Fe Ka line (see Section 3.3). For both sources, two Gaussian functions are 
included in the model, for the narrow Fe K/3 and the broad Fe K emission line, with parameters fixed as in Table [2] Power-law nor- 
malization parameter C is in units of keV -1 cirT 2 s~' at 1 keV; wind column density Nh is given in units of crrT 2 , and the ionisation 
parameter £ in units of erg cm s~'. The two-phase warm absorber parameters are constrained within the range evaluated with the 
double Comptonisation model (Table [3J. The last two columns give the ratio of the chi-squared value to the number of degrees of 
freedom and the associated probability. 
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